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Abstract 

Using the kinctic-theory-basod two- fluid models as a starting point, wc develop 
filtered two-fluid models for a gas-particle flow in the presence of an isothermal, 
first-order, solid-catalyzed reaction of a gaseous species. As a consequence of 
the filtering proc:edure, terms describing the filtered reaction rate and filtered 
reactant dispersion need to be constituted in order to close the filtered species 
balance equation. In this work, a constitutive relation for filtered reaction rate is 
developed by performing fine-grid, two-fluid model simulations of an isothermal, 
solid-catalyzed, first-order reaction in a periodic domain. It is observed that the 
cluster-seale effectiveness factor, defined as the ratio between the reaction rate 
observed in a fine-grid simulation to that observed in a coarse-grid simulation, 
can be substantially smaller than unity, and it manifests an inverted bell shape 
dependence on filtered particle volume fraction in all simulation cases. More- 
over, the magnitude of the deviation in the cluster- scale effectiveness factor 
from unity is a strong function of the m,eso-scale Thielc modulus and dimen- 
sionless filter size. Thus coarse-grid simulations of a reacting gas-particle flow 
will overestimate the reaction rate if the cluster-scale effectiveness factor is not 
accounted for. 
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1. Introduction 



Gas-particle flows are commonly encountered in the chemical and energy 
conversion industries in the form of fluidized beds, risers, and other pneumatic 
conveying units. The flow behaviors of gas-particle systems are commonly an- 
alyzed using continuum models that treat the particle and fluid phases as in- 



terpenetrating continua (Gidaspow 1994 Fan and Zhu 1998 Jackson 2000). 
Finite- volume simulations of these continuum, or 'two-fluid' models, reveal that 
the scale of the gas-particle flow structures that are observed is a strong function 
of the grid size used in the simulation, with grid size independent results being 
obtained when the ratio between the grid size and particle size is 0(10) ( [Agrawal 



et al. 


2001 


Andrews IV et al. 


2005 


Igcietal. 2008 



tide diameters are not computationally affordable when performing two-fluid 
model simulations of large fluidized bed reactors that are often tens of meters 
tall and many meters in cross-section. As a result, coarse-grid simulations of 
large fluidized beds are common practice. However, these coarse-grid simula- 
tions effectively neglect the presence of fine scale gas-particle flow structures 
that are manifested by the two-fluid model. To enable accurate coarse-grid sim- 
ulation of the two-fluid model equations one must account for the effect of the 



fine-scale gas-particle flow structures in the coarse-grid simulation ( Andrews IV 



et al. 


2005 


Igci et al. 


2008 



velopment of a filtered two-fluid model approach for non-reacting gas-particle 
flows where effective closures for the meso-scale fluid-particle drag force, particle 
phase stress, and particle phase viscosity were extracted from fine-grid simula- 



tions of the two-fluid model equations (Igci et al. 2008 Igci and Sundaresan 



2011a|b[ |Igci et al.[ |2011[ [Parmentier et al.[ |2011[ ). While filtered two- fluid 



models for non-reacting monodisperse gas-particle flows have been developed, 
verified against fine-grid simulations, and validated against experimental obser- 
vations (Igci and Sundaresan 2011a|b Igci et al. 2011), the extension of these 
models to reacting and polydisperse gas-particle flows remains an open problem. 
Fine-grid two-fluid model simulations of reacting gas-particle flow have been 
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shown to yield quantitative agreement with experimental observations of con- 
version in solid-catalyzed ozone decomposition processes when grid-independent 



solutions are obtained ( Syamlal and O'Brien 2003 ). However, due to the compu- 
tational expense of performing grid-independent simulations, coarse numerical 
resolution is often used to simulate reacting gas-particle flows in the continuum 
model framework. By not accounting for the fine-scale structure these coarse- 
grid simulations have been shown to over-predict the conversion of ozone in 



a bubbling fiuidized bed (Zimmermann and Taghipour 2005). In addition, a 
few studies have also demonstrated that small-scale gas-particle flow structures 
that form in devices like fluidized beds are responsible for the wide variability 



in mass transfer models that is present within the fluidization literature (Dong 
et al.[|2008a"|bl|Kashyap and Gidaspow||2010[[20TI| ). 



The primary objective of this work is to demonstrate the need for filtered 
models for reacting gas-particle flows by performing fine-grid simulations of a 
first-order, isothermal, solid catalyzed reaction in a periodic domain and filtering 
the results. It will be shown that the presence of particle clustering in fine-grid 
simulations leads to an effective reaction rate that is substantially smaller than 
what would be predicted via coarse-grid simulation. We define the ratio of the 
reaction rate in the fine-grid simulation to that in the coarse-grid simulation 
as the cluster-scale effectiveness factor. It will be shown that this cluster-scale 
effectiveness factor is a strong function of dimensionless filter size A = A\g\/v'^, 
and other model parameters. Here, A is the filter size, g is the gravitational 
acceleration vector, and vt is the terminal settling velocity of an isolated particle. 

Finally, it is shown that grid resolutions finer than those used to deduce fil- 
tered models for non-reacting gas-particle fiows are necessary in order to obtain 
grid-independent filtered closures for reacting gas-particle flows. This observa- 



tion is supported by a recent work by Cloete et al. (2011 ). Due to this grid size 
dependence we propose a filtered reaction rate model based upon extrapolated 
effectiveness factor data obtained from fine-grid simulations. 
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2. Microscopic two-fluid model equations 

The two-fluid model consists of balance equations for gas- and particle-phase 
mass and momentum with an additional transport equation governing evolution 
of the particle-phase granular temperature. The evolution equations that con- 
stitute this Eulerian framework are presented in Table 1, in addition to the 
constitutive relations that close the balance equations. The constitutive rela- 
tions given in Table 1 were chosen to be consistent with earlier studies in filtered 



model development within our research group (Agrawal et al. 2001 Andrews 



|1V et al.[ [20051 feci et all [20081 llgcTand SundaresaSl |2011b[ |Igci et al.[ [20TTj 



Igci and Sundaresan 2011a). Here we note that while we have restricted our 
attention to a specific set of constitutive relations, it has been shown that the 
clustering and bubbling phenomena manifested by the two-fluid model frame- 



work are robust to changes in constitutive relations (Glasser et al. 19981. How- 
ever, we do expect that the constitutive relations comprising the microscopic 
two-fluid models will have quantitative effects on the filtered models that arise 
from fine-grid simulations. Therefore, while we have employed a certain set of 
constitutive relations in this study, filtered models can be developed from fine- 
grid, two-fluid model simulations using a different set of constitutive relations 
as well. 

The effective diffusion coefficient D* and effective gas-phase viscosity /i* 
appearmg m eqs. ^ and ^ will differ fr om the corresponding molecular prop- 
erties as they also account for the enhanced scalar and momentum transport 
that occurs due to additional pseudo-turbulent transport occurring as a result 
of microscopic interactions between individual particles and the fluid. However, 
two-fluid model simulations have revealed that the solutions manifested by the 



two-fluid model are insensitive to the value of /i* (Agrawal et al. 2001 ). There- 
fore, we set ^* equal to the molecular viscosity jjig. In addition, the effective 
diffusion coefhcient D* is kept on the order of the bulk molecular diffusivity D, 
which is generally 0(10^^) rr? j sec for many different gas species. 

Here we investigate the role of clustering on a model, isothermal, solid- 
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catalyzed, first-order chemical reaction 



B Ri^ -kcS(t)PgXgA (1) 

where (j) is the particle volume fraction, and XgA mass fraction of species 

A in the gas phase. For first order reaction kinetics, the effective reaction rate 
constant /ccff based on the bulk concentration Xgi is related to the intrinsic 
reaction rate constant k via 

^ 77»fc ^ 1 f 1 _ 

l + $2^,/Sz $ ^tanh(3$) (3$)/' 

where r]i is the intra- particle effectiveness in the absence of mass transport limi- 



tations, <i> = ka/Dj is the Thiele modulus, a is the volume to surface area ra- 



tio, and Bi is the Biot number for mass transport given as Bi = kmu/Dj (Raw! 



ings and Ekerdt 2009). Here, km is the convective mass transport coefficient, 
a is the particle radius, and Di is the intra-particle diffusivity. It is impor- 
tant to note that the effective reaction rate constants used in two-fluid model 
simulations incorporate the combined effect of mass transport resistance and 
intra-particle effects as well. In this study all gas-solid flow computations are 
performed for FCC catalyst particles fluidized by air — the physical properties 
for solid and fluid phases are given in Table 2. 



3. Filtered tviro-fluid model equations 

The filtered two-fluid models are obtained by performing a spatial average 
of the microscopic two-fluid model equations. The consequence of the filtering 
approach is that the fine-scale gas-particle flow structure that occurs on a length 
scale smaller than the filter size is accounted for through residual terms that 
must be constituted from theoretical considerations or fine-grid two-fluid model 
simulation results. In the case of non-reacting monodisperse gas-particle flows, 
filtered models with accompanying constitutive relations for residual terms have 
been shown to yield quantitatively similar macroscopic behaviors to those ob- 
served in fine-grid simulations of the same fiow problem, with a dramatic savings 
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Evolution Equations 

d{Ps(t>) 

dt 



+ V • {ps(l>v) = (3) 



+ V • {pgegu) = (4) 

dipgGgXgi) 



dt 

d (ps^v) 



+ V • {PgegXg.u) = -V • (£>*egVXsJ + (5) 



+ V • {Ps(t>Vv) = -0V • CTg - V • o-s + /d + Ps(t>g (6) 



3 " 



2 

Closures 



h V • ipsfpTv) 



+ V • (Pgeguu) = -CgV a-g- fn+ PgCgg (7) 
-V q-CTs -.Vv- Jvis - Jcoii + Tsiip (8) 



Gas phase stress tensor (<Tg) 

ag=PgI-p;(vu + {Vuf -^{V-u)I^ (9) 
Solid phase stress tensor (cTs) 

= (p,<^(l + 4<^5o)r-W6(V •«))/- 2m.S, '7=^^^, (10) 
1 2f>QiJ.(fPgQ bp.dVwT . 
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Fluid-particle drag force (/r>) 



S=l(Vv + {Vvf)-^-{V-v)I. (13) 



fn=Hu-v), p=^Reg4>plef -'^ ^ e,p,d|n .;| ^^^^ 

, 18 (Re-^ + 0.15i?e-°-3i3) ije„ < 1000 

Cd = { ^ " ' ' " (15) 

0.33 Rcg > 1000 
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Diffusive flux of granular energy (q) 



q = -AVT (16) 

A = £ ( (l + f ^05o) (l + f r?^ (4r? - 3) <^5o) + |^ (41 - 33r?) ^V^5o' 

(17) 



X* _ Aj _ 75psdVTTT 

5{Ps4>) 9oT 

Dissipation of granular energy via inelastic collisions ( Jcou) 

Jeon = ^ry(l-r,)^<?oT3/2 (19) 
Dissipation of granular energy via viscous action of the fluid phase (Jvis) 

Jyis = ^ -Rdiss (20) 

| + — </.ln(,^)+11.26</. (1 - 5.1.^ + 16.57.^2 _ 21.77<i)^)-^go In(O.Ol) 



(21) 



Production of granular energy through slip between phases (rgUp) 



^ slip ~ _ _ ryf^ ^ , -3 K n: , K nj. ^'^^> 



god^PgV^ 1 + 3.5^0 + 5.9(/) 

l+3A/^+(135/64)<^ln((/.)+17.14<^ 



< 0.40 



_ I 1+0.6810-8.4802+8.1603 ^23) 



Table 1: Microscropic two-fluid model equations for an isothermal, reacting gas-particle flow. 
Here, ps and pg are particle and fluid phase densities respectively; (p is the particle volume 
fraction; €g is the gas void fraction; v and u are particle and fluid phase velocities; Ri is 
the rate of production of species i; D* is the effective diffusion coefficient; //* and fj,s are the 
effective gas phase viscosity and solid phase viscosity, respectively; 50 is the radial distribution 
function at contact; Cp is the coefficient of restitution; A is the conductivity of pseudothermal 
energy; and T is the granular temperature 



7 



d 


7.5 X IQ-^ m 


Ps 


1500 kg/m^ 


Pg 


1.3 kg/m^ 


Pa 


1.8x10^5 kg/ {m s) 




0.9 



Table 2: List of physical parameters for solid and fluid phases. 



in computational time (Igci and Sundaresan 2011b I. Here we follow, the filter- 



ing procedure given in Igci et al. ( 2008 ) . 

The particle volume fraction and gas-species mass fractions obtained from 
fine-grid simulations can be represented as (j){y, t), where y and t represent the 
location and time variables, respectively. The filtered particle volume fraction 
{(j)) {xq, t) is given as 

{^)ixo,t)= f G{xo,y)cj>{y,t)dy (24) 

where G{xQ,y) is a weight function, and Xq is the spatial location of the filter 
center. In this work (•) will be used to indicate filtered variables. We require 
that jy G{xo, y)dy = 1, and in this study we use a top hat filter for G{xo,y). 
Filtered gas-species mass fraction (xg;) {xq, t), particle phase velocity (v) {xq, t), 
and fluid phase velocity (u) {xo,t) are given as 



(eg) {xo,t) (xgj {xo,t) 

((/)) {xo,t) (v) {xo,t) = 
(eg) {xo,t) (u) {xo,t) = 



G{xo,y)egiy,t)xg,{y,t)dy (25) 
G{xo,y)c^{y,t)v{y,t)dy (26) 



Gixo,y)eg{y,t)u{y,t)dy. (27) 

Upon filtering the microscopic species balance equation for component A we 
obtain the following filtered species balance equation 
diPa (Cff) (Xqa)) 



dt 



+ V • (pg {egXg^u)) = -V-{D (egVxg^)) - W<, ('/'Xg J • (28) 



In eq. (28) we invoke the definition of filtered variables given in eqs. (25l~(27). 
Due to the fact that the reaction we are considering in this work is isothermal 
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and produces no volume change, the filtered momentum balance equations re- 



main unchanged from those derived in the work of Igci et al. ( 2008 ) , and as 



such, they will not be presented here for the sake of brevity. In cq. (281, there 



are several terms that appear as filtered products of microscopic variables that 
must be constituted to solve the filtered two-fluid model equations. To facili- 
tate filtered model simulations we must decompose the products of microscopic 
variables into mean and fluctuating parts. Here, we define fluctuating variables 
as follows 

Xg^ (y, t) = (y, t) ~ ixg^) {y, t) <) = t) - (u) (y, t) (29) 

0'(y,t) = (/.(y,t)~(0)(y,i) (30) 

where x'g^ i , and 0' represent species A mass fraction, gas velocity, and volume 
fraction fluctuations, respectively. Inserting the decomposition of mean and 



fluctuating parts into eq. (281, we obtain the following flltered species balance 
equation 



d {Pg {eg} iXg.)) ^ ^ . ^^^^ ^^^^^ («))=- V . {pg (e.x'gy) + D (e,Vx,J) 



dt 



kcSPg (((/') (XgA) + (Xg^)) 



(31) 



The flrst term appearing on the right hand side of eq. (31 1 can be modeled 



as a dispersive term given by the following constitutive equation 

^cffV ixg,) - {pg (egXgy) + D (e.Vx.J) , (32) 
where D^s is a dispersion coefficient. Dispersion coefficients defined in this 



way have been presented in the research literature (Loezos and Sundaresan 



20021. The final term on the right hand side of eq. (31) represents a filtered 



reaction rate that must also be constituted in terms of filtered variables alone. A 
straightforward method for constituting this reaction rate is to define a cluster- 
scale effectiveness factor rj/^ which is given as 

filtered reaction rate {'t'Xgi) 



VA = 



homogeneous reaction rate {(p) [xgi 



(33) 
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However, in this paper we present filtered models and associated constitutive re- 
lations for the non-locally corrected cluster-scale effectiveness factor 77^ defined 
as 

where TO2 = A^/12. The approach of removing non-local effects in the determi- 



nation of filtered quantities was recently advanced by Parmentier et al. (2011), 
and we apply the same procedure here for determining rj'^. The non-locally 
corrected effectiveness factor is constructed by removing the dominant gradient 
terms that contribute the evaluation of the filtered product {(pXgi), which is 
done to allow the construction of filtered constitutive relations in terms of local 
filtered variables alone. We present a derivation of the non-locally corrected 
effectiveness factor in [Appendix A| It will be shown below that both 77A and 
77^ depend on volume fraction in a qualitatively similar way, but differ quanti- 
tatively. All data presented in future sections and resulting filtered models will 
pertain to rj'^, thus necessitating one to track the gradients in filtered species 
mass fraction and volume fraction to enable the use of the resulting filtered 
models in coarse-grid simulations of reacting gas-particle flow. 

Dimensional analysis of the parameters governing gas-particle flow suggests 
that the cluster-scale effectiveness factor is a function of five independent di- 
mensionless quantities, which are 

n, = y^ = $ n2 = ^ n3 = i^ = 5~c (35) 

V D vt PgD 

n^^^ n^^^^A. (36) 

Pg -"t 

One can also readily define additional dimensionless parameters like Reg = 
PgVt/iig and the coefficient of restitution Cp using dimensional analysis. How- 



ever, earlier works of Igci et al. (2008) have shown that filtered quantities do 



not display any significant dependence on Rcg or Cp, and as such, they were not 
included in the dimensional analysis presented here. Moreover, we found that 
H2 and n4 had a negligible effect on the cluster-scale effectiveness factor, while 
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the meso- scale Thick modulus $, Schmidt number iSc, and the dimensionless 
filter size A substantially alter the cluster-scale effectiveness factor. Therefore, 
all work presented below will only interrogate the effect of Sc, and A on 
77^. Here the term meso-scale Thiele modulus is used to distinguish it from the 
Thiele modulus based on intra-particle difFusivity Di given in eq. ([2]). 

4. Numerical implementation 

When performing periodic domain simulations of reacting gas-particle flows, 
one is faced with the following dilemma: due to the imposition of periodic 
boundary conditions, the concentration of any reactant undergoing an irre- 
versible reaction within the periodic domain decay to zero. In order to facilitate 
gas-particle flow simulations in periodic domains, an alternate simulation strat- 
egy must be developed. For the special case of a first-order reaction, it can 
readily be shown that one can track the evolution of — Xg. / (Xgi) rather than 
Xgi itself, and relate the filtered value of Hi directly to the cluster-scale effec- 
tiveness factor. The benefit of tracking the evolution of Ki lies in the fact that 
it has a non-zero statistically steady value, even though Xgi will decay to zero. 

Consider taking the species balance equation given by eq. ([5| with the re- 
action rate expression given by eq. ([T]) and averaging it over the entire periodic 
domain, resulting in the following equation 

P9 {^9) %^ ^-V. f kom<l>X9,dV, (37) 



V 



dt V 

where V is the volume of the periodic domain. If we now define the variable 
= Xgi I {Xgi ) E^nd plug Kj into the species balance equation given by eq. ([5| 
we obtain the following evolution equation 



d{Pg(-gHi) 



V • {pgEgKiu) = - V • {DegS/ni) " hsPg^iKi 

KsPgcjiKidV. 



{^g)V 

Due to the fact that the reaction rate expressions in this problem are first order, 
we are able to track the evolution of without considering the time progression 
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of Xg. or {xgt) ■ The presence of non-zero source and sink terms on the right hand 



side of eq. ( 38 1 force statistically steady value of Ki to be non-zero regardless of 
the reaction rate constant used. Therefore, in all work presented in following 
sections for the evolution of Ki is solved, and the cluster-scale effectiveness factor 
is redefined in terms of 

All simulation results presented in this work were generated using the Multi- 
phase Flow with Interface eXchanges (MFIX) software that relies on a variable 
time step, staggered grid, finite- volume method for the solution of the two-fluid 



model equations (Syamlal 1998). The iterative solution to the two- fluid model 
equations is obtained using the Semi-Implicit Method with Pressure Linked 
Equations, or the SIMPLE algorithm. Due to the strong coupling between solid 
and fluid phases that arises due to the fluid-particle drag force, the Partial Elim- 
ination Algorithm of Spalding was used to effectively decouple the solution of 



the solid and fluid phase balance equations (Spalding 1980). In addition, a sec- 
ond order Superbee discretization is employed for the convective terms that are 
present in the conservation equations for continuity, momentum and granular 
energy transport to limit the effects of numerical diffusion. 

5. Preliminary Study 

To illustrate the behavior of the cluster-scale effectiveness factor as a func- 
tion of {(f)) , simulation results are presented that were obtained using grid sizes 
that are sixteen times as large as the particle diameter. This grid resolution 



was chosen to follow the simulations employed in the study of Igci et al. (2008). 
Filtered variables were obtained from these fine-grid simulations by moving fil- 
ters of different sizes throughout the periodic domain. Due to the statistical 
homogeneity of periodic domains we are able to collect thousands of samples 
and average them regardless of spatial position. Moreover, by filtering these 
fine-grid simulation results for various domain-averaged volume fractions the 
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filter-averaged volume fraction dependence of the cluster-scale effectiveness fac- 
tor can be ascertained. This filtering procedure follows that outhned in the 



presented in Figure [Tj It is clear from Figure [T] that both rj^ and rj'^ are strong 
functions of (0) retaining an inverted bell shape that approaches unity in the 
limit of small and large particle volume fraction. However, there are noticeable 
quantitative differences between rj^ and 77^ in Figure^with the maximum devi- 
ation being near the minimum in both curves. The minimum value of 77A ~ 0.33, 
while the minimum value of 77^ w 0.53. Therefore, at this resolution, the non- 
local correction to the cluster-scale effectiveness factor contributes as much as 
37% to the value of the cluster-scale effectiveness factor near the minimum in 
the curves given in Figure [T] While the quantitative differences between rjA 
and 77^ are functions of A and other model parameters, this example provides 
a scale of the difference between the two effectiveness factors. In this work we 
have decided to model 77^ because it is directly related the small scale fluctua- 
tions in (0) and (xg;), while r/A is dependent on non-local variations in {(f>) and 
(Xsi) (see Appendix). Physical intuition for the decrease in the cluster-scale 
effectiveness factor from unity can be obtained by observing the characteristic 
clustering patterns that are observed at different domain-averaged volume frac- 
tions, which are superimposed in Figure [l] At low volume fractions there are 
a few small isolated clusters throughout the periodic domain, and as such, the 
cluster-scale effectiveness factor will begin to deviate from unity. As the vol- 
ume fraction increases the frequency of these clusters increases thus making the 
effective contacting between gas and particle phases poor. Near the minimum 
value in the 77^ curve, clusters begin to span the periodic domain, and the solid 
phase changes from the dispersed to the continuous phase. As volume fraction 
continues to increase the gas-particle flow becomes more homogeneous, and due 
to this homogeneity the cluster-scale effectiveness factor increases toward unity 
in the limit of high particle volume fraction. 

In Figure [2] (a) the cluster-scale effectiveness factor is shown as a function 
of {(p) for four different $ values at fixed values of the Sc and A. The depth of 




dependence of 77^ and 77^ is 
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0.2 0.4 0.6 



Figure 1: A characteristic plot of the 77^ versus (</<) is shown with snapshots of the particle 
volume fraction field different domain-averaged volume fractions superimposed. Here, <I> = 
0.23, Sc = 1.00, and A = 3.14. 
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the inverted bell shape curve in {(j)) is an increasing function of $. The values 
of $ given in Figure [2] (a) are substantially smaller than unity, but there is a 
marked change in rj'^ from unity. This may seem peculiar, but the small values 
of $ arise as a result of the fact that the length scale used to determine $ 
is the particle diameter d. A more fitting length scale is that associated with 
a cluster. However, the size of a cluster emerges as a result of fine-grid two- 
fluid model simulations, and cannot be considered as an input parameter. It 
is for this reason that we choose the particle diameter as the relevant length 
scale for $ rather than the length scale of a cluster. In Figure [2] (b) the {(f) 
dependence of the cluster-scale effectiveness factor is given for four different 
values of the dimensionless filter size A at fixed $ and Sc. The departure of 
the effectiveness factor from unity is also an increasing function of A. However, 
Figure [2] (c) shows that the change in the cluster-scale effectiveness factor from 
unity is a decreasing function of Sc when keeping $ and A constant. While 
one might expect that increasing the Sc value should increase the departure in 
rj'^ from unity, we observe a decrease in this departure due to the fact that we 
are demanding $ to remain constant in Figure [2] (c). In these simulations the 
variation of Sc is achieved by varying the value of D. In order to maintain a 
constant $ the effective rate constant k^s must be changed. Therefore, varying 
Sc at constant <i> requires the variation in the effective rate constant. It is this 
coupled variation that brings about the somewhat surprising the dependence of 
rj'^ on Sc. 

6. Grid Resolution Effect 

When deducing a filtered model from a fine-grid simulation, it is necessary 
to ensure that filtered statistics are independent of grid size 5. To that end, the 
dependence of B = 1 — r;^ on (</>) is presented for four different grid resolutions 
in Figure [s] (a) at fixed 4, A, and Sc. The dependence of if^ on grid size 5 
is substantial, and only begins to saturate when the grid size is around 4 — 8 
particle diameters. Figure 2 (b) shows the variation of the dimensionless filtered 
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Figure 2: The cluster eflfectiveness factor r;^ as a function of filtered particle volume fraction 
(</>). (a) rj'^ versus {(p) for different values of the meso-scale Thiele modulus $ for fixed values 

of the meso-scale Schmidt number Sc and dimensionless filter size A; (b) r;^ versus {(p) for 
difi'erent values of A for fixed values of Sc and "I"; and (c) rj'^ versus {rj>) for different values of 
Sc for fixed values of A and 
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A = 1.57 



1.57 




(a) 



O S/d = 32 
S/d = 16 
> A/d = 8 
□ <5/d = 4 




0.2 0.3 



(b) 



Figure 3: Grid size <5 dependence of filtered quantities, (a) B = 1 — r)'^ is plotted versus (0) for 
several different grid resolutions, (b) The dimensionless filtered drag coefficient {f}e.,d) is plot- 
ted versus (</<) for several different grid resolutions to show the degree to which hydrodynamic 
variables are fully resolved. Here, <I> = 0.23 and Sc = f .00. 



drag coefRcient (/Je.d) as a function of (0) for four different grid resolutions. 
Here the dimensionless filtered drag coefficient is defined consistent with the 



earlier work of Igci et al. ( 2008 ) 



vt ({fD,)-{4>'dp'Jdvy 



Ps\g\ 



(40) 



where the subscript y indicates that the drag coefficient is inferred from the 
drag force and pressure fiuctuation terms directed parallel to gravity. At grid 
sizes of sixteen particle diameters {Pe.d) seems to exhibit grid independence for 
{(p) < 0.30, with some grid dependence emerging at higher volume fractions. 
Figures 3 (a) and (b) illustrate the sensitive grid dependence of the cluster-scale 
effectiveness factor when compared to other filtered quantities like the fiuid- 
particle drag coefficient. Consequently, even finer grid resolutions are required 
when simulating reacting gas-particle flow when compared to non-reacting sys- 
tems. This observation further supports the need for the development of filtered 
models for accurate coarse-grid simulation of reacting systems. 

Due to the sensitivity of 77^ with respect to grid size, we seek to develop 
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a model for the cluster-scale effectiveness factor extrapolated to the limit of 
infinite grid resolution. In Figures |4] (a) and (c), plots of B/ Bmax versus {4>) /(p* 
are given for the finest two grid resolutions presented in this study. Here, B^ax 
is the maximum value of B, and (jf is the filtered volume fraction at which 
the maximum value of B occurs. From inspection of Figures |4] (a) and (c) it 
is clear that the grid dependence of the volume fraction variation in B can be 
removed by rescaling B by B^^ax and (0) by (jf for all (0) < cf)* . In addition, 
it is demonstrated in Figures |4] (b) and (d) that the grid dependence of the 
volume fraction variation in B can be removed for all (0) > 0* by plotting 
B/Bjnax versus {{(p) — cj)*) / {(f>c — (f>*)- Here (j)c is the volume fraction at which 
the value of B reaches zero. Due to the grid-independent behavior observed 
in Figures |4] (a) — (d) a filtered model utilizing a piecewise description of the 
variation in B / B^ax with particle volume fraction about cf)* is developed, while 
extrapolating the values of B^^ax: 4'* ^ and (p^ to infinite resolution. 

In Figure [5] (a) and (b) the grid size dependence of B„iax is presented for 
two different values of $ and a variety of different filter sizes. The values of 
Bmax are clearly saturating as the value of 5 approaches zero. In order to pro- 
vide an accurate value of B^ax to use in our filtered reaction rate model, the 
value of Bmax at infinite resolution is determined via a Richardson Extrapola- 



tion (Roache 1998). The values of 0* and 0c were observed to vary linearly 
with grid resolution, and a linear extrapolation was performed to ascertain the 
infinitely resolved estimates of </>* and 4>c- 



7. Extrapolated results and filtered model 

In the previous section, grid independence of the cluster-scale effectiveness 
factor was observed by plotting B/Bmax against scaled volume fraction coor- 
dinates that differ depending on whether (0) is greater than or less than <j>* . 
Motivated by this observation, we seek to model the volume fraction variation 
of B/Bmax via a piecewise function about </>*. In Figures [6] (a) — (c) plots of 
B/Bmax versus ((/)) /(j)* are presented for diflferent values of <l>, A, and Sc. The 
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Figure 4: Here we demonstrate the grid size independence of the cluster-scale eflfectiveness 
factor curves when appropriately scaled, and the volume fraction dependence is considered 
separately to left and right of (f)* . In (a) and (c) we plot B = 1 — rj'^ scaled by its maximum 
value Bmax as a function of {rj>) /(/>* for the finest two grid resolutions in this study. In (b) 
and (d) we plot B/Bmax as a function of {{(p) — </>*) / {4>c — <t>*)- 
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Figure 5: We present B^ax plotted versus S/d for several different values of A. (a) "I = 0.11 
and Sc = 1.00; (b) I- = 0.16 and Sc = 1.00. 



shape of the curve of B /Bmax is clearly altered by changes in $ and Sc, with 
only mild changes in the volume fraction dependence of B/Bmax at different 
grid values of A. In Figure|6](d) the variation in B/B„iax with (0) /ip* is shown 
to collapse by plotting all data with the same value of ^^A/Sc together. 

Figures [7] (a) — (c), show the variation of B/B„iax with {{(/)) — <p*) / (0c — 0*) 
for different values of $, A, and Sc. The Sc variation in the dependence of 
B/Bjnax on {{(p) — (f>*) / {(pc — 0*) is clear, while the $ and A dependence is sub- 
stantially weaker. In Figure[7](d) we illustrate that the dependence of B/B^ax 
on ((0) — (j)*) I (jpc — 0*) for different <I>, A, and Sc values can be collapsed onto 
a single curve provided the value of $'^A/S'c is kept constant. From Figures [g] 
(d) and[7](d) it is clear that the relevant scaling parameter governing the shape 
of the dependence of B/B,nax on (0) on both sides of (j)* is dictated by a single 
parameter given as $^A/S'c. The following functional forms are used to describe 
the {(f)) dependence of B/B,nax- 

B ^( l-{l-{<t>) /<j>*f ^ 

Bmax \l + ai{l- {(p) /(t)*f J 



(0) < r 



(41) 
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Figure 6: Dependence of B/B.max on {(f)) /0* for varfous values (a) (b) A, and (c)Sc. In 
(d) we present B/Bmax versus {(j>) /<p* where 9^K/Sc is constant. 
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Here, Oi and a2 are least-squares fit parameters that depend only on <i&^A/S'c, 
see Figures |8] (a) and (b). Both model parameters are given by smooth functions 
of <i>^ A/ Sc presented in Tablejs] It should be noted that the observed fluctuation 
in a2 about the model curve in Figure |8](b) can be attributed to the uncertainty 
in determining (pc from the simulation results for 77^. 

In Figure [9] (a) Bmax is presented as a function of A for a variety of different 
4> and Sc values. The variation of Bmax with $, A, and Sc is evident by 
inspection of Figurejo] (a). However, by replotting all data as a function of <i>^A 
alone one can collapse all B^ax values onto a single master curve in Figure |9] 
(b). A curve fit is presented for B„iax in terms of $^A in Table 3. 

Linear extrapolation of 0* to infinite resolution reveals that (jj* varies be- 
tween 0.34 <(/)*< 0.50, with no discernable trend in $, A, or Sc. This variation 
in (j)* arises as a result of the fact that the cluster-scale effectiveness factor is 
a weak function of (0) in the region around (j)* and thus determining the value 
of the (j)* is subject to error. Since no discernable trend in (jj* was observed, we 
recommend the use of cj)* — 0.42 because it represents the ensemble average of 
the different 0* values obtained. Moreover, since B only varies slightly in the 
vicinity of B^ax, small errors in <j)* will not influence the quantitative behavior 
our filtered model substantially. In addition, the extrapolated value of (pc was 
found to vary between 0.58 < 4>c < 0.62 with no systematic dependence on $, 
A, or Sc. As a result of this fluctuation we choose (pc = 0.59 because this is 
consistent with earlier work in our group suggesting that filtered model correc- 
tions for gas-particle hydrodynamics become negligible at volume fractions of 



0.59 and higher (Igci and Sundaresan 2011a 



The cluster-scale effectiveness model developed in this work predicts that 
effective reaction rates observed in coarse-grid simulations be larger than those 
predicted in fine grid simulations if the effects of fine scale structure are not 
accounted for via the cluster-scale effectiveness factor. Therefore, in order to 
accurately perform continuum model simulations of reacting gas-particle flows 
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Figure 7: Dependence of B/Bmax on ((</)} — </)*) / (c/ic — </>*) for various values (a) <E>, (b) A, 
and (c)Sc. In (d) we present B/Bmax versus ((</>) — </>*)/ (0c — </>*) where 9^K/Sc is constant. 
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Figure 8: Fit parameters describing the dependence of B/Bmax with {<!>). (a) Gives the 
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on coarse spatial grids, filtered models for effective reactions rates are neces- 
sary. Without these corrections, coarse-grid continuum model simulations will 
consistently over-estimate the conversion observed in a solid catalyzed reaction 



(for example see Zimmermann and Taghipour (2005)). Therefore, in order to 
perform accurate coarse-grid numerical simulations of a first-order, isothermal, 
solid-catalyzed gas-phase reaction we suggest the use of the following filtered 
species balance equation 

dip, ^ ^ . ^^^^ ^^^^^ («» = - V • (i^effV - KspWa {<!>) (x.J 

(43) 

where the constitutive relation for rj'^ is given in Table |3j and a constitutive 
relation for the effective dispersion coefficient D^k can be inferred from the work 
of Loezos and Sundaresan (2002 ). In section[5j it was shown in that the non-local 
correction to the effectiveness factor contribute much as 37% to the observed 
value of 77a- We have observed that this contribution decreases as a function of 
increasing grid resolution and filter size. Therefore, in the limit of large filter 
sizes we expect that the contribution of the non-local term {(f))^ • V (xg.) 



in eq. (43) will be weak, and can thus be neglected. However, for intermediate 
filter sizes we recommend the inclusion of this non-local correction. 

While in this work the evolution of species mass fraction was solved in the gas 
phase alone, one could envision a model for a solid-catalyzed, gas phase reaction 
where the species mass fraction in both gas and solid phases is tracked separately 
and coupled through a mass transport term between particle and fluid-phases. 
In a model framework of this type, the effective interphase mass transfer rates 
will be decreased due to the presence of clustering in the gas-particle flow, and 
any reduction in reactant conversion that arises can be attributed to decreased 
mass transfer efficiency. We expect that such an observed decrease in mass 
transfer coefficient will be on the order of the cluster-scale effectiveness factor 
developed in this work. As an example, consider a mass transfer operation 
taking place between FCC particles and air, where the effective mass trans- 
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4.4V12A 

7^ (44) 

1 + 4.4V$2A 



1 -^max 



2 7fi 

«i - . . (46) 

exp n7.2$2A/5'cj 

02 = 265$2a/S'c + 10.4 (47) 



Table 3: Filtered model for the cluster-scale effectiveness factor valid for a first order, solid 
catalyzed, isothermal reaction. 



fer coefficient is calculated to be /c,„ w 0(0. 1) m/sec using a model developed 
by Gunn (19781 for fixed beds. Using A;,„ « 0.1 m/sec one can determine the 
characteristic rate constant for mass transfer k' = kmSp/Vp, where Sp and Vp 
are the surface area and volume of an FCC particle, respectively. For an FCC 
particle k' « 8000 sec~^, and using the duster-scale effectiveness factor model 
developed in eq. (44) for a filter size of 2 cm we predict that 77^ = 0.07. Us- 
ing this minimum value of the cluster-scale effectiveness factor we predict that 



fcgff — rj'^k' ~ 530 sec 



a nearly 20 fold decrease in the effective mass trans- 



fer coefficient! Indeed decreased effective mass transfer coefficients on this order 
have been found in Energy Minimization Multi-Scale Model (EMMS) simula- 



tions of mass transport processes in gas-particle flows ( Dong et al. , 2008a ) . 



8. Conclusions 

The need for the development of filtered two-fluid models for reacting gas- 
particle flows is demonstrated by considering a model first-order, isothermal 
solid-catalyzed gas-phase reaction. It is shown that constitutive relations for the 
filtered reaction rate and filtered species dispersion must be postulated to close 
the filtered species balance equation. Due to the fact that the model gas-phase 
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reaction in this work is isothermal and produces no volume change, only filtered 
species balance equations must be developed, without the need to alter the 
existing filtered gas and solid momentum balance equations given in the earlier 



works of Igci et al. (20081 and Igci and Sundaresan (2011a) for non-reacting, 
monodisperse gas-particle flows. Using flne-grid continuum model simulations 
of reacting gas-particle flows we extract the cluster-scale effectiveness factor, 
defined as the ratio of the fine-grid reaction rate to the reaction rate observed in a 
coarse grid simulation. The cluster-scale effectiveness factor is observed to retain 
an inverted bell shaped dependence on volume fraction approaching unity in 
both the low and high volume fraction limits. At intermediate volume fractions 
a decrease in the cluster-scale effectiveness factor from unity is observed with 
the magnitude of this reduction being a strong function of meso-scale Thiele 
modulus <f>, and dimensionless filter size A with only a weak dependence on 
Schmidt number Sc. 

Due to the sensitivity of the cluster-scale effectiveness factor with grid size, 
we determine an asymptotic form for the cluster-scale effectiveness factor re- 
lying on a Richardson extrapolation of the minimum in the cluster scale effec- 
tiveness factor. The extrapolated values of the cluster-scale effectiveness factor 
are shown to collapse when plotted as a function of <1>^A, with the volume 
fraction dependence of the cluster-scale effectiveness factor collapsing for fixed 
^^A/Sc. Tabled presents a curve fit of our collapsed results for use in coarse- 
grid simulations of reacting gas-particle flows with first-order reaction kinetics. 
The reduction in effective reaction rates observed in this study can be used 
to rationalize the over-prediction in reactant conversion that was seen in the 



two- fluid model simulations of Zimmermann and Taghipour ( 2005 1 . Finally, 
we note that the filtered model developed in this work is limited to first-order 
reaction kinetics, and was restricted to two-dimensional periodic domain simu- 
lations. Future work should extend such analyses to other reaction kinetics, and 
three-dimensional bounded domains. However, we expect that the characteristic 
volume fraction, Thiele modulus, filter size, and Schmidt number dependence 
of 77^ will be qualitatively similar to the results presented here. 
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9. Notation 



Cd - fluid-particle drag coeflacient a - particle radius 

ai - fit parameter for plot of B/Bmax versus {(j)) /</>* 

a2 - fit parameter for plot of B/Bmax versus ((</)) — (j)*) / {(pc — 0*) 

B-B = l-r,'^ 

Bi - Biot number for mass transport Bi = kmCi/D 

Bmax - maximum value of B d - particle diameter 

D* - effective diffusivity of reacting species in continuum model 

D - molecular diffusivity of reacting species 

Dj - intra-particle species diffusivity 

Des - meso-scale dispersion coefficient 

Bp - coefficient of restitution of the particle phase 

fo - fluid-particle drag force experienced by the fluid 

g - gravitational acceleration vector 

go - radial distribution function at contact 

G{xo,y) - weight function for filtering 

■/coll - coUisional dissipation of granular energy 

Jvis - viscous dissipation of granular energy 

k - intrinsic reaction rate constant 

k' - effective rate constant for mass transfer 

fceflF - effective reaction rate constant 

fcgjf - meso-scale effective rate constant for mass transfer 

km - convective mass transport coefficient 

m2-m2 = A2/12 

Pg - gas-phase pressure 

q - granular energy conduction vector 

Ri - rate of production of species i 

Reg - Reynolds number for the gas-phase 

Sp - Surface area of a particle 

S - rate of deformation tensor 
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Sc - meso-scale Schmidt number Sc = Hg/ (pgD) 
T - granular temperature 
u - gas velocity 

V - particle velocity 

Vt - terminal settling velocity of an isolated particle 

V - volume of a periodic domain 
Vp - Volume of a particle 

xq - generic spatial position vector 
xq - spatial position of filter center 
y - generic spatial position vector 

Greek letters: 

j5 - fluid-particle friction coeSicient 

Tsiip - production of granular energy through interphase slip 

5 - dimensional grid size 

A - dimensional filter size 

A - dimensionless filter size h. = /S.\g\/v1 

eg - gas volume fraction 

?? - r? = (1 + ep) /2 

r}i - intraparticle eff'ectiveness factor 

77a - cluster-scale effectiveness factor 

T}'^ - non-locally corrected cluster-scale effectiveness factor 

Kj - ratio of Xgi to its domain- averaged value 

A - conductivity of granular energy 

/Xb - bulk viscosity of particle phase 

fig - molecular gas-phase shear viscosity 

li*g - effective gas-phase shear viscosity 

lis - shear viscosity of particle phase 

Pg - gas density 

Ps - particle density 

ag - gas-phase stress tensor 
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Us - particle-phase stress tensor 

- particle volume fraction 

$ - meso-scale Thiele modulus $ ~ ^JTzflJJdFjD 

$ - Thiele modulus 4 = y/k(P/Di 

Xgi - mass fraction of gas species i 
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Appendix A. EfTectiveness factor with non-local contributions re- 
moved 

In section [3] the cluster scale effectiveness factor is defined in two ways given 



by eqs. (|33|— (|34|). Here, the effectiveness factor given in eq. (34) is derived 



relying on the definition of the effectiveness factor given in eq. (33). Let the 
particle phase volume fraction and species mass fraction Xgi at £^ny point be 
represented as follows 

(l){y,t) = {cj)){y,t) + cf>'{y,t) Xg.{v.t) = {Xg.){y.t) + ^gAy.t). (A.l) 

where y is a spatial variable associated with filtered and microscopic variables, 
respectively, neither of which are located at the filter center. Assuming the 
filtered variables can be given by smooth functions of space the filtered value of 
4> and Xgi can be approximated via the following Taylor series 

(0) (y, t) = (0) {xo,t) + {y-xo)-\7 + 0{{y - xo) ■ {y - xo)) (A.2) 

iXa.) iy,t) = iXg^) {xo,t) + {y-xo)-V {xg,)\(^^„^t)+Oi{y'Xo)-{y-xo)). (A.3) 

The Taylor expansions given in eqs. (A.2 ) and ( A.3[ ) can be plugged into eq. ( A.l[ ) 
to yield expressions for cj) and Xgi that depend only on the location of the filter 
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center Xq and the microscopic spatial variable y. Utilizing eqs. (A.l) — (A.3) the 
filtered product of (f) and Xgi can be expressed as 



(0Xg.) ixo,t) ^ {(f)) {xo, t) (xg,) (xo, t) 

+ m2V(0)^- V(x3j|(xo,t) 

^'{y,t)x'g^{y,t)G{xo,y)dV 



(A.4) 



where TO2 — A^/12. Removing the gradient term that appears on the right hand 



side of eq. (A.4) provides a method to define filtered constitutive relations for 



77^ that depend on local filtered quantities alone, without considering gradients 
in local filtered variables. This approach of removing non-local effects in the 
process of filtered model development has recently been advanced when con- 
sidering the development of filtered fluid-particle drag models for monodisperse 



gas-particle flows ( Parmentier et al. 2011 ), and this formulation is extended here 



to reacting gas-particle flows. Combining the definition of 77A and 77^ given in 



eqs. (33 1 and (34) with eq. (A.4) one can arrive at the following two expressions 



(r/A - 1) (0) ixo,t) (xg.) (xo, t) 



(t>'{y,t)x'{y,t)G{xo,y)dV 



+ m2V • V(X5J (Xo,t)- 



(A.5) 



(77k-l)(0)(a;o,t)(XgJ(a^o,i)= / q^'{y,t)x'{y,t)Gixo,y)dV. (A.6) 



Jv 

Therefore, 77^ provides a direct measure of the product of local fluctuations in 
{(p) and (Xgi)i while the value of ryA is influenced by gradients in local filtered 
quantities. In order to interrogate the effect of local fluctuations alone, we 
choose to model rj'^ in this study. 
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